Hermitian and non-Hermitian formulations of the time evolution of quantum decay 
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This work discusses Hermitian and non-Hermitian formulations for the time evolution of quantum 
decay, that involve respectively, continuum wave functions and resonant states, to show that they 
lead to an identical description for a large class of well behaved potentials. Our approach is based 
on the analytical properties of the outgoing Green's function to the problem in the complex wave 
number plane. 
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I. INTRODUCTION 

The theoretical description of quantum decay refers to 
the time evolution \ipt) = exp(— iHt/h)\ipo) of an initial 
state |-0o) m a system characterized by a Hamiltonian H . 
In some decay problems it is not convenient to separate 
the Hamiltonian into a part with stationary states and 
a part which is responsible for the decay, usually treated 
to some order of perturbation, but rather to consider the 
full Hamiltonian H to the system. This is usually the 
case when the decay originates by tunneling through a 
classically forbidden region. As is well known, follow- 
ing the work by Khalfin [1], if the energy spectra E of 
the system is bounded by below, i.e., E 6 (0, oo), the 
exponential decay law cannot hold at long times. At 
short times there is also a departure from the exponen- 
tial decaying behavior which is related, however, to the 
existence of the energy moments of the Hamiltonian H 
[2-4] . These both type of behaviors have been confirmed 
experimentally in recent times [5, 6]. 

The motivation for this work goes back to the early 
times of quantum mechanics. In 1928, Gamow intro- 
duced the notion of resonant state to describe the time 
evolution of decay of a particles in radioactive nuclei [7] . 
In order to describe the above process, Gamow consid- 
ered solutions to the Schrodingcr equation which at large 
distances consist only of purely outgoing waves. This is 
physically appealing because it yields an outward flux for 
the decaying particle outside the interaction region. He 
realized, however, that the absence of incoming waves 
in the solution at large distances leads to complex en- 
ergy eigenvalues. This was in a way satisfactory because 
it led to the interpretation of the imaginary part of the 
energy as the inverse of the lifetime r in the exponen- 
tial decay law exp(— t/r), and thus it provided a the- 
oretical framework for the understanding of the expo- 
nential decay law in quantum mechanics. In fact, his 
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approach constituted one of the first successful appli- 
cations of quantum mechanics involving tunneling phe- 
nomena and also one of the first theoretical treatments 
of open quantum systems. However, since the ampli- 
tude of resonant states increases exponentially with dis- 
tance, the usual rules of normalization, orthogonality and 
completeness do not apply. Nonetheless, over the years 
in spite of these apparent drawbacks, a consistent theo- 
retical framework involving resonant states evolved over 
the years [8, 9]. The approach involving resonant states 
represents a non-Hermitian formulation that lies strictly 
outside the usual Hermitian framework of quantum me- 
chanics. In fact, there is a traditional view that considers 
the non-Hermitian description of quantum decay involv- 
ing complex energy eigenvalues as an approximate, phe- 
nomenological description that cannot be fundamental 
because it violates the requirement of unitarity [10]. The 
formulation of decay involving continuum wave functions 
requires of numerical integration over the wave number 
for each value of the time [11, 12]. It turns out, however, 
that the approaches mentioned above may lead to re- 
sults for the time evolution of decay that are numerically 
indistinguishable from each other [13]. So the question 
naturally arises of how these very different formulations 
of decay lead to the same numerical results. 

The aim of this work is to investigate the genesis of the 
above two formulations, Hermitian and non-Hermitian, 
for the decay process using the analytical properties of 
the outgoing Green's function to the problem in the com- 
plex k plane. We intend to answer or at least to throw 
some light in understanding the question posed at the 
end of the above paragraph. 

The organization of the paper is as follows. In Sec. 
II some general properties of the time-dependent wave 
solution involving the outgoing Green's function to the 
problem are briefly discussed. Subsection II A refers to 
the continuum wave solutions and provides a derivation 
for time-dependent solution in terms of continuum states. 
Subsection II B yields a discussion on the formalism of 
resonant states and yields an expression for the time- 
dependent solution in terms of resonant states. In sec- 
tion III, we illustrate equivalence between the basis of 
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continuum and resonant states for the J-shell potential. 
Finally, Sec. IV gives the concluding remarks. 



II. TIME-DEPENDENT SOLUTION AND THE 
OUTGOING GREEN'S FUNCTION 

We shall consider a very simple, yet no trivial, descrip- 
tion of the time evolution of decay of a particle that is 
confined initially in a real spherical potential of arbitrary 
shape in three dimensions. Without loss of generality we 
restrict the discussion to s waves and it is worth men- 
tioning that the description holds also on the half-line 
in one dimension. It is also assumed that the interac- 
tion potential V(r) vanishes after a finite distance, i.e. 
V(r) = for r > a. However, as discussed below, the re- 
sults obtained hold also for potentials that go faster than 
exponentials at very large distances. The units employed 
here are h = 2m = 1. 

Let us therefore write the time-dependent Schrodinger 
equation as 



Im k 



,d_ 

di 



H 



*(r,t) = 0, 



(1) 



where the Hamiltonian H = —d 2 /dr 2 + V(r). 
The solution to Eq. (1) may be written as 



*(r,t) 



g(r,r';t)*(r',0)dr', t>0 (2) 



where g(r,r';t) stands for the retarded time-dependent 
Green's function to the problem. This function describes 
the time evolution of the system for t > and has van- 
ishing value for t < [14] . Clearly, in order to obtain the 
time-dependent solution one must know g(r, r ; t). A con- 
venient form to determine this quantity is by expressing 
it in terms of the outgoing Green's function to the prob- 
lem G + (r,r';k). For t > 0, both quantities are related 
by using the Laplace transform method [15] 



<?(r,r';*) = -l_ f G +{r,r';k) 

27T« J C 



-ik t 



2kdk, t > 0, 



(3) 

where Co, before taking the radius of the semicircle Cr 
up to infinity, corresponds to the contour in the complex 
k plane shown in Fig. 1. One may then deform the 
contour as shown in that figure. The exponential factor 
in the integrand of (3) guarantees that the contour Cr 
along the second quadrant of the k plane vanishes as 
the radius of the semicircle goes to infinity. As a result 
g{r, r';t) becomes 



g(r, r';t) 



i 

2^ 



G+{r, r';k)e 



ik t 



2kdk, t > 0. 



(4) 

In deriving the above expression we have assumed also 
for simplicity, since decay refers to a process where the 
particle tunnels out into the continuum, that the poten- 
tial does not hold bound states. 




Re k 



FIG. 1. Contour Co and its deformation along the complex k 
plane used to derive Eq. (4). 



Equation (4) is the central quantity to study the dy- 
namics of decay and it is the starting point to analyze 
the genesis of the descriptions involving continuum wave 
functions and resonant states. One sees that it depends 
on the outgoing Green's function and hence it is conve- 
nient to refer first to some of its properties. The function 
G + (r, r'; fc) obeys the equation 

[k 2 -H]G + (r,r';k) = S(r-r'), (5) 

and satisfies the outgoing boundary conditions, 

G+(0,r';k) = 0, G'+(a,r';k) = ikG+(a,r';k), (6) 

where the prime denotes here, and thereafter, the deriva- 
tive with respect to r. The outgoing Green's function 
may be written in terms of the so called regular, 4>{k, r), 
and irregular, f±{k, r), solutions to the Schrodinger equa- 
tion and the Jost function as [14] 



G+(r, r'-k) 



H k , r <)f+( k > r >) 
J + (fc) 



(7) 



where r < and r > stand respectively for the smaller and 
larger of r and r' . We shall discuss briefly some properties 
°f f± an( i J± an d then return to discuss Eq. (7). The 
regular solution (f>(k, r) satisfies the equation 



[k 2 — H](f)(k, r) = 0, 



(8) 



with boundary conditions <f>(k, 0) = and </>'(fc,0) = 1. 
The irregular solutions f±(k,r) satisfy 



[k 2 - H]f ± (k,r) =0, 



with boundary conditions 



f±(k,r) = e^ mi , r > a. 
Finally, the Jost functions are defined as 

J±(k) - [f±<t>' - f ± 4>]r=a. 



(9) 



(10) 



(11) 
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Also, J±(k) = f±{k, 0). The regular and irregular solu- 
tions satisfy, respectively, integral equations that are of 
the Volterra type and hence may be solved by iteration 
[14]. This leads to certain conditions on the behavior of 
the interaction potential V(r) as a function of the dis- 
tance r. These arc that the first and second moments of 
V(r) must be finite, which imply, respectively, that the 
potential is less singular near the origin that r~ 2 and goes 
faster than r -3 as r — > oo. In addition, it is required that 
J dr r\V (r)\e nr < oo, with r\ a real quantity. In partic- 
ular, if the potential decreases faster than exponential 
at infinity, /+ and /_ are entire functions of k in the 
whole complex k plane. This is the case considered in 
this work, where the interaction vanishes beyond a dis- 
tance, but holds also for other type of potentials, as for 
example, potentials with Gaussian tails. 

Some relevant properties of <p and f± for real val- 
ues of k are [14]: 4>(k,r) = (f>(—k,r) = (f)*(k,r) and 
f—(k,r) = f1(k,r). Notice that the first relationship 
in the above expressions implies that <f>(k, r) is real and 
an even function of k. The irregular solutions /+ and /_ 
are linearly independent functions and it may be shown 
that <f> may be written as a linear combination of them 

<Kk, = 2^ t J ~ ( fc )/+ ( fc > r ) - J + ( fc )/- ( fc > r )l . ( 12 ) 
which along region r > a, in view of (10), reads 

<j>(k,r) = ^ [J-(k)e lkr - J+(k)e~ ikr ] 
The Jost functions may also be expressed as [14] 



(13) 



J±(k) = l + k~ 1 dr sin krV(r)f±{k,r), (14) 
Jo 

and for real values of k one has J-(k) = J+(k). 

A consequence of the above considerations is that the 
outgoing Green's function G + (r, r';k) is single valued 
and analytical in the whole complex k plane except at 
an infinite number of poles that correspond to the zeros 
of the Jost function J+(fc). For potentials that vanish 
after a distance these poles are in general simple and we 
shall assume that this is the case here. One finds, in gen- 
eral, a finite number of them may seat on the positive and 
negative imaginary k axis, corresponding respectively to 
bound and antibound states, and thatban infinite num- 
ber are located on the lower half of the k plane where, 
due to time-reversal considerations, are distributed sym- 
metrically with respect to the imaginary k axis, corre- 
sponding, as discussed below, to resonant states. Thus 
for a pole n n = a n — i/3 n on the fourth quadrant of the 
k plane, there corresponds a pole «_„ = — «* that seats 
on the third quadrant. 



A. Time-dependent solution in terms of continuum 
wave functions 

Here we obtain an expression for the time-dependent 
solution given by Eq. (2) as an expansion in terms of 



continuum wave functions. The discussion, of course, 
involves values of k that are real. It starts by noticing 
that Eq. (4) may be written as 



(r,r';t) = — / [G+(r, r'; k)-G+(r, r'; -k)]e~ lk l 2kdk 

27T Jn 

(15) 

Then using Eqs. (7) and (12) one may write the inte- 
grand to Eq. (15) as 



2- 



[G + (r 7 r'-k) - G+{r,r';-k)]2k = i: + (k, r)^ + *(k, r'), 

(16) 



where the continuum wave functions arc defined as 



7T J+(fc) 

Substitution of (16) into (15) gives 



7T Jl(k) 



(17) 



g(r,r';t)= V + (fc,r)V + ' {k,r')e- lk 1 dk, t > 0, 
Jo 

(18) 

and substitution of (18) into (2) allows to write the time- 
dependent solution as 

/>oo 

*(r,t)= / C(k)i; + (k,r)e- lk2t dk, (19) 
Jo 

where the expansion coefficient C{k) is given by 

C(k)=[ i/j + *(k,r')^(r\0)dr'. (20) 



Equation (19) yields the time evolution of the time- 
dependent solution ^(r, t) as an expansion in terms of 
the continuum wave functions to the problem. Notice 
that by taking the limit as t —> in Eq. (18) yields the 
closure relationship 



4>+{k,r)iP + *{k,r')dk = S(r -r'), t>0, (21) 



which shows that the set of continuum wave functions is 
complete. 

The continuum wave functions are solutions to the 
Schrodiger equation of the problem 

[k 2 - H]^+(k,r) =0, (22) 

and satisfy the boundary conditions 

ip + (k,Q) = 



V + (fc,r) 



2 i 
n2 



-ikr 



ikrl 



r> a, (23) 



where S(k) is the S-matrix of the problem. A comparison 
of the second term on the right-hand side of (23) and (13) 
allows also to write the S-matrix as 

J_(k) 



(24) 



J+(k) 

and similarly, a comparison between (23) and (12) leads 
to Eq. (17), that relates the continuum wave functions 
with the regular solutions and the Jost functions. 
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B. Time-dependent solution in terms of resonant 
states 

Resonant states are defined as the solutions to the 
Schrodingcr equation 

[«* - H]u n (r) = (25) 

obeying the boundary conditions, 

u„(0) = 0, u' n ( a ) = iK n u n (a). (26) 

Notice that the second of the above conditions means 
that for r > a, u n (r) = D n exp(iK„r), and hence, as 
first discussed by Gamow [7] , it involves complex energy 
eigenvalues k^ = E n = £ n — iT n /2, with K n = a n — i0 n , 
where £ n = c? n - f3 n and T n = Aa n /3 n . 

The modern approach to resonant states is based on 
the analytical properties of the outgoing Green's func- 
tion on the complex k plane. The aim is to obtain an 
expansion of G + (r, r'; k) in terms of its poles. 

Let us therefore consider the expression [16] 



2m J c k' — k 



(27) 



where C is a large closed contour of radius L in the k' 
plane about the origin, which excludes all the poles K n 
and the value k' = k, namely, C = Cr + + ^2 n c n . 
Choosing Cr in the clockwise direction, and Ck and the 
contours c n in the counterclokwise direction, it follows 
using Cauchy's theorem, that J = 0, and hence one may 
write 



2niJ 



E 



Cn 



G + (k>) 
k'-k 



dk' 



G+{k') 
k'-k 



dk' 



G+(k'] 



dk' = 0. 



(28) 



It turns out that as the radius L — ¥ oo, which increases 
the number of poles inside the contour Cr up to infinity, 
the function G + (r, r';k), that appears in the integral over 
the circle Cr, diverges unless both r and r' arc smaller 
than the potential radius a, or one of them has the value a 
and the other remains smaller than a [17, 18]. We denote 
the above conditions, that guarantees convergence of the 
resonant sum in (28), by the notation (r, r')t < a. One 
may then use the theorem of residues to evaluate the re- 
maining terms in (28). This requires to know the residues 
p n at the poles K n of the outgoing Green's function. They 
follow by adapting to the k plane the derivation given in 
Ref. [15], namely, 



Pn(r, r') 



u n (r)u n (r') 



2n n {/ Q a ul(r)dr + iul(a)/2K n } 



(29) 



which provides the normalization condition for resonant 
states 



u^(r)dr + i 



= l. 



(30) 



Notice that for bound states, where Kb = ijb, (30) re- 
duces itself to the usual expression. A orthogonality con- 
dition for resonant states follows using Green's theorem 
for solutions u n and u rn of (25) and its corresponding 
boundary conditions, to obtain 



u n (r)u m (r)dr + t^^^ = 0. (31) 

&n "T Km 



Hence, using (29) in view of (30) gives the purely discrete 
expansion 



G + {r,r';k)= ]T 



71— — OO 



Un{r)u{r' n ) 



(r, r') f < a. (32) 



Substitution of (32) into (5) leads to the following rela- 
tionships 

- Yl u n (r)u n (r') = S(r - r'), (r, r') f < a, (33) 



which stands for a closure relation, and the sum rule 

g u n {r)u n (r') = Qj r/)f ^ a (34) 



Notice that l/[2« n (fc- K n )] = l/(2fc[l/(fc- k„) + 1/k„]), 
and hence (32), in view of (34), may be written also as, 

G+(r,r';*) = ^ £ (r,r'y<a. 

n— — oo 

(35) 

Substituting (35) into (5) yields again (33) and the new 
sum rule 

oo 

J2 u n (r)u n (r')K n = 0, (r, r') f < a. (36) 



n— — oo 



Substitution of (35) into (4) provides the resonant expan- 
sion for the time-dependent Green's function along the 
internal interaction region, namely, 



g (r, r'; t) = ^ u n (r)u n (r')M(y°), (r, r') f < a, 



(37) 



where M(y°) stands for the Moshinsky function [19] 

i r°° e -ik 2 t i 
M(y° n ) = - / -. dk = -w(iy° n ), (38) 

Zn J-oo K - K n I 

with y° = —e^^^Knt 1 / 2 . The function w(iy°) stands for 
the Faddeyeva function [20] for which well known com- 
puting algorithms have been developed [21]. 

One may also derive an expression for G + (r, r'; k) along 
the external interaction region r > a by noticing that for 
r' < a and r > a, G + (r,r';k) = G + (a,r';k)e ik ^ r - a \ 
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Substitution of this expression into (4) and expanding 
G + (a,r';k), using (35), gives 

oo 

g(r,r';t) = ^ u n (r')u n (a)M(y n ), r' < a, r > a, 

(39) 



n— — oo 



where M(y n ) stands for the Moshinsky function 



M(y n 



i 

2^ 



^ik(r—a) p —ik t 



-dk 



(40) 



with y n = e- i7r / 4 (l/4i) 1 / 2 [(r - a) - 2n n t}. 

Substitution of (37) and (40) into (2) allows to write 
the resonant expansion for the time-dependent solution 
as 



*M) = E 



where 



C 



C n u n {r)M{y° n ), r<a 
C n u„(a)M(y n ), r>a; 



u n (r)^(r,0)dr. 



(41) 



(42) 



Notice in (41) that at r = a, y n = 

It is worth mentioning that using the properties of the 
Faddeyeva function [20], one may derive an explicit an- 
alytical expression of the time-dependent function along 
the internal region of the potential at long times [8] 



i£„* e -r„t/2 



n=l 



J ^ C n u n (r) 1 1 

I,Im \S~) ^ 



r < a, 



(43) 



where n = l/(47r) 1 ' 2 . Notice that the sum rule (35) 
has canceled out exactly the leading t~ x / 2 term in the 
asymptotic long-time behavior of the Faddeyeva function 
[8,20]. 

The time-dependent decaying solution given by the 
first term in (41) satisfies the continuity equation 
i(d/dt)\^(r,t)\ 2 + (d/dr)J(r,t) = 0, with J(r,t) the cur- 
rent density, which by integration along the internal in- 
teraction region yields the relationship, 

r - 2a Ma)|2 (44) 

In deriving the above expression it is required the deriva- 
tive of the Faddeyeva function [20] and Eqs. (31) and 
(36). 



III. MODEL 

In order to illustrate the Hcrmitian and non-Hcrmitian 
formulations discussed above, we consider a <5-shcll po- 
tential, of radius a and intensity A for s waves, namely, 



V(r) = XS(r - a). 



(45) 



o.o 
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FIG. 2. Behavior of the probability density in lifetime units 
at the potential boundary r — a for the <5-shell potential with 
parameters as indicated in the text. 
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FIG. 3. Behavior of the probability density in lifetime units 
at the potential boundary r = 25a for the <5-shell potential 
with parameters as indicated in the text. 



This model was initially considered by Winter [22], and 
since then, by many authors. The reason being that its 
mathematical simplicity does not avoid that it describes 
correctly the main physical features of the time evolution 
of decay. As initial state we choose the infinite box state, 



\E-(r,0) 



2\ 1/2 (nr\ 
- sin I — . 
a I \ a 1 



(46) 



It follows then, using (23) that the continuum wave func- 
tions, corresponding to the Hcrmitian formulation, read 



?/> + (fc, r) 




sin(fcr)/ J_)_(fc), r < a 



i/2) [e 



— ikr 



S(k)e 



' kr 1 



(47) 



r > a, 



where the S-matrix is given by (24). The Jost function 
J+(k) of the problem may be obtained immediately using 
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(14), namely, J+(k) = 2ik+X(cxp(2ika— 1), and we recall 
that J-(k) — J±(k). One then may obtain expressions 
for the expansion coefficients C(fc), given by (20), and 
then, by using (47), either along the internal or external 
regions of the potential, of the time-dependent solution 
^(r, t) given by (19). To obtain the time-dependent so- 
lution requires of numerical integration along a span of 
values of k for each value of time t. 

On the other hand, for the non-Hcrmitian formulation, 
the resonant states of the problem read, 

( A n sin(K„r), r <a 
u n (r) = { (48) 
[ B n e lK " r , r>a. 

From the continuity of the above solutions and the dis- 
continuity of its derivatives with respect to r (due to 
the (^-function interaction) at the boundary value r = a, 
it follows that the k„'s may be obtained by solving the 
equation, 2in n + A(e 2lK " a — 1) = 0, which corresponds 
to the zeros of the Jost function J + (k) . For A ^> 1 one 
may write the approximate analytical solutions to the 
above equation as, k„ (7r/a)(l — 1/Aa)— i (l/a)(7r/Aa) 2 . 
Using the above expression for n n as the initial value 
in the Newton-Rapshon method, i.e., = K r n — 

F(<)/F(<), with F = [dF/dk] k = Kn yield the solutions 
n n with the desired degree of approximation according 
to the number of iterations. 

The normalization coefficients of resonant states may 
be evaluated by substitution of Eq. (48), for r < a, into 
Eq. (30), to obtain A n = [2A/(Aa + cxp(-2m„a))] 1 / 2 . 
Similarly, using Eqs. (46), the first expression in (48) 
and the above expression, yields analytical expressions 
for the coefficient C n . Hence, for given values of the po- 
tential parameters, A and a, one may then calculate the 



[1] L. A. Khalfin, Sov. Phys.-JETP 6, 1053 (1958). 

[2] L. A. Khalfin, JETP Lett. 8, 65 (1968). 

[3] J. G. Muga, G. W. Wei, and R. F. Snider, Europhys. 

Lett. 35, 247 (1996). 
[4] G. Garcfa-Calderon, V. Riquer, and R. Romo, J. Phys: 

Math. Gen. 34, 4155 (2001). 
[5] C. Rothe, S. I. Hintschich, and A. P. Monkman, Phys. 

Rev. Lett. 96, 163601 (2006). 
[6] S. R. Wilkinson, C. F. Bharucha, M. C. Fischer, K. W. 

Madison, P. R. Morrow, Q. Niu, B. Sundaram, and M. G. 

Raizen, Nature 387, 575 (1997). 
[7] G. Gamow, Z. Phys. 51, 204 (1928). 
[8] G. Garcfa-Calderon, Adv. Quant. Chem. 60, 407 (2010). 
[9] G. Garcfa-Calderon, AIP Conference Proceedings 1334, 

84 (2011). 

[10] C. M. Bender, Rep. Prog. Phys. 70, 947 (2007). 
[11] D. A. Dicus, W. W. Repko, R. F. Schwitters, and T. M. 
Tinsley, Phys. Rev. A 65, 032116 (2002). 



set of complex poles {n n }, resonant states {u n (r)} and 
coefficients {C„}, to evaluate the time dependent solu- 
tions (41). Notice that all these quantities are calculated 
only once, which saves considerable computation time. 

Figures 2 and 3 provide plots of the In \ ijj(r, t) as a func- 
tion of time in lifetime units, r = 1/Ti, for parameters of 
the (5-shell potential, A = 12.0 and a = 1, respectively, for 
r/a = l and r/a = 25. One observes in both cases that 
the Hermitian and non-Hermitian formulations arc indis- 
tinguishable from each other. Also plotted, in both fig- 
ures, is the asymptotic long-time contribution that goes 
as t~ 3 . 



IV. CONCLUDING REMARKS 

It is worth stressing that the exact Hermitian and 
non-Hermitian formulations for the description of decay 
discussed here, corresponding to coherent (elastic) pro- 
cesses, have as common origin the analytical properties 
of the outgoing Green's function to the problem and that 
its equivalence holds for potentials having finite first and 
second moments and tails that go faster than exponen- 
tial at long distances or that vanish after a distance. This 
might allow to consider artificial quantum systems as ul- 
tracold atoms [6] and resonant tunneling structures [23]. 
Since decay by emission of particles corresponds to an 
open quantum system, it is not surprising that the notion 
of unitarity usually employed for closed systems does not 
apply. In this context, Eq. (44) is very relevant because 
it implies flux conservation as time evolves. Finally, it is 
also worth emphasizing that the resonant formalism pro- 
vides expressions, as (43), which exhibits explicitly the 
exponential and noncxponcntial contributions to decay in 
contrast with the 'black box' description (19) involving 
continuum wave functions. 
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